Quasi-continuum approximation to the Nonlinear Schrodinger equation with 

Long-range dispersions 
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The long-wavelength, weak-dispersion limit of the discrete nonlinear Schrodinger equation with 
long-range dispersion is analytically considered. This continuum approximation is carried out irre- 
spective of the dispersion range and hence can be assumed exact in the weak dispersion regime. For 
nonlinear Schrodinger equations showing finite dispersion extents, the long-range parameter is still 
a relevant control parameter allowing to tune the dispersion from short-range to long-range regimes 
with respect to the dispersion extent. The long-range Kac-Baker potential becomes unappropriate 
in this context owing to an "edge anomaly" consisting of vanishing maximum dispersion frequency 
and group velocity (and in turn soliton width) in the "Debye" limit. An improved Kac-Baker po- 
tential is then considered which gives rise to a non-zero maximum frequency, and allows for soliton 
excitations with finite widths in the nonlinear Schrodinger system subjected to the long-range but 
finite-extent dispersion. 

PACS numbers: 03.50.-z, 42.65.Tg, 45.20.-d 



Considerable efforts have been devoted to understand- 
ing physics of one-dimensional(lD) systems in which dy- 
namical properties are dominated by the competition 
between nonlinearity and dispersion. The Nonlinear 
Schrodinger equation(NLSE) is one of most investigated 
nonlinear equations in these contexts, and widespread 
applications can be found in various fields of condensed 
matter physics. Very recently, a new path was opened in 
the interest to this equation toward materials with long- 
range(LR) dispersions These LR dispersions are 

thought to occur in two different kinds namely, intrachain 
dispersions involving long-range interactions among par- 
ticles of a single discrete chain 0] , and interchain disper- 
sions related to the couplings between several ID chains 
with short-range interactions ■ Each of these two kinds 
of LR dispersion deserves interest since both lead to two 
distinct but real physical contexts. Magnetic spin chains 
and Davydov's molecular chains 0, 1^ are two exam- 
ples of natural systems in which the first kind of LR 
dispersion arises. As for the second, we strongly sus- 
pect it may provide an excellent way to model systems 
made of several weakly coupled ID chains. Indeed, if 
the coupling between the ID chains is not too strong 
the model can stand as a good approximation of a quasi- 
one-dimensional((51Z)) material. To this viewpoint, we 
expect the second kind of LR dispersion to find direct 
applications in nonlinear optics or Bosc-Einstcin conden- 
sate(BEC) systems [llj, |12| where soliton compression 
phenomena H, IE likely to occur, e.g. resulting 
from periodic arrangements of dielectric elements as in 
photonic crystals 8, 9, 10] or from the configuration of 
an optically coupled array of BEC lattices. 
The present work deals with the first kind of LR disper- 
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sion. The discrete NLSE in this cas writes 0: 

L 

-ilpnt + ^ Jrn-n{i^ni " V'n) + 5 I V'n | ^ "^n = 0. (1) 

where is the complex wavefunction, 17 is a nonlinear 
coupling(hereafter assumed positive) and Jm-n is the 
potential creating the LR dispersion. Equation ^ has 
been discussed by few authors assuming exponential 
and power-law LR potentials. These works focus on soli- 
ton solutions as well as on their stability and suggest 
rich physical properties related to the LR dispersion. In 
particular, authors of ref. H combined lattice discrete- 
ness and LR dispersion and showed that different soli- 
ton regimes may be stabilized from the interplay of these 
two factors. In connection with the variety of stability 
regimes, the authors suggested a switching mechanism by 
which the transition between different bistable localized 
states could be monitored by tuning parameters of the 
LR potential. In biophysical systems where inter molecu- 
lar interactions often display great sensitivity to the flex- 
ibility of molecular backbones, such a switching mecha- 
nism is likely to be due to the " contraction- relaxation" 
features of bond danglings. In this view, provided the 
bond length is commensurate with the wavelength of in- 
tra molecular excitations, the discrete NLSE with LR 
dispersion(LRNLSE) is indeed a good approximation for 
nonlinear molecular vibrations and nonlinear molecular 
excitons can well be understood as highly discrete and lo- 
calized solutions of this equation with relatively narrow 
shapes(compared with the lattice constant). However, if 
the bond length is short and the molecular masses large, 
narrow excitons are less probable since long-wavelength 
processes will tend to dominate. This last situation is 
the commonly studied one in magnetic spin chains and 
Davydov models in the short-range dispersion regime. 
The associate NLSE admits single pulse and dark soli- 
ton solutions in the continuum limit. Strictly speak- 
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ing, the continuum short-range NLSE(SRNLSE) is the 
long- wavelength and weak-dispersion limit of the discrete 
SRNLSE 13]. 

The goal of the present study is to point out that 
an equivalent limit is also accessible for the discrete 
LRNLSE irrespective of the range of lattice dispersion. 
For this purpose, we will follow a Fourier transform 
method consisting to first construct the dispersion func- 
tion of the LR system, and next extract the continuum 
NLSE from a second-order expansion of the dispersion 
function. As this expansion does not involves some con- 
straint on the range of lattice dispersion, the continuum 
LRNLSE derived from our long-wavelength and weak- 
dispersion approximation is indeed an exact equivalent 
of the celebrated continuum SRNLSE. Our interest to 
the problem was motivated after noting that the disper- 
sion curves of the NLSE with exponential LR potentials 
always displayed the linear shape at small wavectors, and 
that only their slopes were affected by the variation of the 
LR parameter. In this context, no physical or mathemat- 
ical constraint prevents us from linearizing the dispersion 
relation of the LRNLSE so far as the long-wavelength 
limit is concerned. Actually, the variation of the slope 
of the dispersion curve with varying dispersion range is 
the signature of a varying sound speed. That the en- 
tire effect of the LR dispersion is confined on the sound 
speed agrees quite well with the spirit of weak-dispersion 
approximations. Disorder, whcih is another interesting 
source of residual dispersion, often causes the same ef- 
fect. We would also like to draw attention on similarity 
between the current problem and that discussed by Kono- 
top 01 in terms of the group velocity of a discrete non- 
linear monoatomic lattice with Kac-Baker(KB) IT^ 
LR interaction. 

For illustrative purpose, we will assume two different 
LR potentials one of which is the KB potential i.e.(^ = 
m — n): 

JKB{e)^Jo^^r\'\, 0<r<l, (2) 
2r 

The second LR(MKB) potential is defined as: 

Jmkb{1)^Jo ^ }~\. r\'K 0<r<l (3) 

In both ||2Jl and ||3Jl, Jo is the potential constant, r is 
the LR parameter and L is half of the spatial extent of 
the LR dispersion. Indeed, we suppose the actual dis- 
tance on which LR dispersions effectively spread out is 
not absolutely equal to the total size of the spin or molec- 
ular system. For instance we can think of an infinite-size 
system undergoing finite-extent dispersions as it really 
happens in most cases. Neverthess, to avoid losing track 
of the usual necessary constraints on LR dispersion po- 
tentials it is needful to mathematically probe the validity 
of such possibility. For this purpose we sum Jm-n over 



L obtaining: 

L 

jRBii) = Jo(l-r^), 

L 

J2 JMKBii) - Jo. (4) 

As one sees, the norm of the KB potential is length 
dependent and becomes constant only in the "thermo- 
dynamic" limit. The dependence of the total ampli- 
tude of the KB potential on L, which we term "edge 
anomaly" [l7| . is the main motivation of our interest to 
the MKB potential. Of course, this edge anomaly does 
not shades at all the rich physics behind the KB model 
though seeming to restrict its validity to infinite-length 
materials and systems with infinite dispersion extents. 
Still, making the potential amplitude a function of L as 
in the MKB potential can turn into a relevant advantage 
for NLSEs with finite dispersion extents. In this case the 
potential amplitude can be adjusted to a desired disper- 
sion extent regardless of whether the total length of the 
system is finite or not. To this point, note that Jmkb 
reduces to Jo in the short-range limit(i.e. L — 1) irre- 
spective of the value of r. On the contrary, for Jkb we 
need to set r = 0. Below a finite-length solution of the 
continuum LRNLSE will be derived and will appear to 
be a periodic function, the period of which can be set 
equal to L. Doing so, we confine solitonic excitations 
on the effective spatial interval spread by the dispersion. 
Since the LR paramater r allows monitoring the strength 
of the dispersion regardless of its extent, the theory will 
still be valid whether the system is of finite size or not. 
As we are interested with the nonlinear localized solution 
of eq. 1^, we can write down the wavefunction ipn{t) as 
follow: 

i'nit) = (j)n exp{-ivt). (5) 
Inserting jSJ in |^ we obtain: 

L 

-V(t)n + ^ Jm-n{4>,n - <l>n) + 94>n = 0. (6) 

Since the new wavefunction is only space dependent, © 
is a discrete nonlinear static equation and a spatial 
Fourier transform can be performed. We define the 
Fourier transform(per unit length) of 0^ as: 

4>n = J dqct>qe-'''^'', (7) 

where d is the lattice spacing and q the wavector. Except 
the LR potential, all coefficients in eq. © are constant 
so their Fourier transforms are single inegrals. Remark 
that the nonlinear term of this equation describes an 
isotropic "three- mode" coupling and as consequence the 
associate Fourier tranform will be a simple integration 



FIG. 1: Plot of the dispersion function versus ^(in the re- 
duced unit (//tt) for some values of the long-range parameter 
r. Here the extent of dispersion is taken infinite. 

over the single wavector common to the three coupled 
fields. The most interestmg part of the discrete nonlin- 
ear static equation is thus the second term whose Fourier 
tranform writes: 

oc „ 
^ J™-„('^™ - = - J dq Mr) <j^q e-^«"^ 

L 

Ml) = J{r)J2r'[l-cos{qM)], 

(8) 

where J{r) can assume two distinct forms in connection 
with lO and © i.e.: 

1 - r 

JKB{r) = Jo , 

r 

1 - r 

JMKsir) = Jo—r-, FT- (9) 

The quantity Jq(r) is the function governing the spatial 
dispersion of the LRNLSE. The sum appearing in this 
function is restricted to L. However, it can be extended 
to infinity provided the chain length also is infinite. Since 
this last context is more familiar, let us first examine the 
case L ^ oo. On fig. ^ we plot Jq versus q for some 
values of r. The just mentioned sum is evaluated numer- 
ically for the two LR potentials. According to |(2Jl, both 
potentials reduce to the same expression in the infinite- 
dispersion limit. 

Fig^ shows that Jq{f) remains a sinusoidal function 
for the whole interval of values of r, and only slopes of 
the dispersion curves vary. To check this assertion, a 
numerical expansion of the dispersion function 3q was 
carried out up to the second-order term for weak values 
of q d. The leading and first-order terms vanish but not 
the second-order term. The coefficient of this last term is 
nothing else but the square of the sound speed Co(r) 
which is plotted in fig. [3 as a function of the LR param- 
eter r. 



FIG. 2: Plot of the sound speed versus r for the NLSE with 
infinite dispersion extent. 



An exact analytical expression of Jq{r) for finite L can 
be obtained using the following identity: 

L 

Mq) ^ '^/cosiqM), 

1=1 

r cos{q) — — A^+'^'> [cos((L + l)q) ~ r cos{Lq)\ 
1 — 2r cos{q) + 

(10) 

and is given by: 

Mq)^Jir)[Vr{0)^Vr{q)], (11) 

We can check that this analytical expression is fully con- 
sistent with numerical curves in fig. ^ To see the effects 
of finite L on the dispersion of our system, in figs. |3| 
and 0] we have drawn Hll() versus r placing ourselves at 
the edge of the first Brillouin zone i.e. q = tt, and using 
L = 10, 20, 1000. We see that the maximum frequency 
is sensible to the extent of dispersion for the KB poten- 
tial and not for the MKB potential. However, the two 
sound speeds(figs.ElandEJl fill effects of finite values of L 
but tend to the same limit when L becomes sufficiently 
large(i.e. infinite). In the case of KB potential, there 
occur maxima in Cl(?') at r < 1 and the sound speed 
vanishes as r ^ 1 for finite L. As opposed to this first 
model, Cl (r) is monotonous in the MKB case and tends 
to a finite value as r ^ 1 when L assumes finite values. 

The analytical expression of the sound speed plotted 
in figs El and El was derived from (|11|) by an analytical 
expansion of Jq{r) to the second order for weak values 
of the product q d. Since here too the leading and first- 
order terms vanish, we are only left with the second-order 
term proportional to q^. Strictly, this expansion cor- 
responds to a continuum-limit approximation for long- 
wavelength(n(i x) and weak-dispersion (qd << 1) ex- 
citations. Thus, with regard to the minus sign in front of 
the right-hand side of the first equation in © , an inverse 
Fourier tranform of the q^ term produces the following 



FIG. 3: Dependence of the maximum frequency on the dis- 
persion extent for the LRNLSE with the KB potential. 



FIG. 5: Dependence of the sound speed on the dispersion 
extent for the KB potential. 




continuum nonlinear second-order equation: 



(12) 



The square of the sound speed C\{r) is proportional to 
the second-order derivative of Jq{r) i.e.: 

CUt) = (l/2)a,,r!2(q)|,^o. (13) 

It is useful to note that the discrete LRNLSE itself 
could be treated in a similar manner and would lead to 
the usual form of continuum NLSE. 
For positive values oig and equation H12|l admits finite- 
length soliton solutions given in terms of Jacobi El- 
liptic functions [l^ : 



X 

,cn I — |K 



_2if_ 

2k2 - 1 \g 



K > 1, 

1/2 



il = ^^{2k^ - 1) 



(14) 



cn, which is the Jacobi Elliptic function, is periodic in x 
with a period Lo = AtoK where K is the complete Elliptic 



integral of the first kind. In general, cnoidal waves pro- 
vide adequate representations of the soliton solutions of 
finite-length systems in connection with the finite magni- 
tude of Lo- In this viewpoint, we can decide to set L = Lq 
which amounts to confine all solitonic excitations within 
the interval spread by the LR dispersion. In the limit 
K = 1, H14|) turns to the single-pulse soliton: 



{x) = iposech 



X 

To 



Cl{r) 



(15) 



and Lo — > oo. However, L can be kept finite in this limit. 
According to the dependence of £o on the sound speed, 
we expect the pulse width to assume behaviours of Cl (r) 
with respect to the LR parameter r. On fig. El Coo{r) 
was an infinitely increasing function of r. This means the 
pulse width is quickly diverging as r — > 1 for infinite dis- 
persion extents. Following the effect of an assumption of 
finite values of L on the sound speed, finite-width single- 
pulse solitons are likely to be observed in a NLSE with 
LR dispersion and corresponds, in a physical context. 
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to a system in which the extent of dispersion docs not 
foUow the system size. However, since the sound speed 
always goes to zero as r 1 for the KB potential, this 
model is relatively less suitable in such physical contexts 
as opposed to the apparently advantageous features of 
the MKB potential. 
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